Determining a quantum state by means of a single apparatus. 
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The unknown state p of a quantum system S is determined by letting it interact with an auxihary 
system A, the initial state of which is known. A one-to-one mapping can thus be realized between 
the density matrix p and the probabilities of occurrence of the eigenvalues of a single and factorized 
observable of S-f A, so that p can be determined by repeated measurements using a single apparatus. 
If S and A are spins, it suffices to measure simultaneously their z-components after a controlled 
interaction. The most robust setups are determined in this case, for an initially pure or a completely 
disordered state of A. They involve an Ising or anisotropic Heisenberg coupling and an external field. 



Consider a set of identical quantum systems S, pre- 
pared in some unknown state p. Information on p can be 
gathered by measuring some observable (I; of S. Such a 
measurement has a statistical nature, whether p is pure 
or not: when repeated many times, it provides the prob- 
abilities Pi = trpTTi for the distinct eigenvalues uji of 
an operator 6j, where tt^ are the associated eigenprojec- 
tions. We thus find partial information on p. A question 
then arises which lies at the heart of quantum theory: 
Which observables need to be measured for a complete 
determination of p iS*^? The recognition that certain 
non-commuting observables have to be measured for that 
purpose was used by Bohr to formulate the principle of 
complementarityi*^. Later on the problem of determining 
an unknown state was considered from various perspec- 
tives for continuous^ and discrete systems^ and found ap- 
plications in quantum communicationSii. However, dur- 
ing all these developments it was not questioned whether 
non-commutative measurements are truly needed. 

Here we face the problem of determining from some 
collection of experimental data the whole set of unknown 
matrix elements of the state /5 of a system S. In fact, 
what we call "the state of a system" refers to a density 
matrix describing an ensemble, and the determination of 
the data will require repeated experiments. However, we 
will show that it is sufficient for our purposes to pro- 
duce these data by means of a single apparatus which 
measures commuting observables only. The key of the 
method consists in coupling the system S to an auxiliary 
system A, the state of which we know. We will show 
that the full p can then be deduced from simultaneous 
measurements of two obviously commuting observables, 
to pertaining to S and 6 pertaining to A, respectively. 
Repeated measurements of lj and 6 yield the joint prob- 
abilities of occurrence for all pairs of eigenvalues of lu and 
o; this will be sufficient to determine the whole density 
matrix p of S, provided the initial state of A and the 
joint evolution of S+A are known. We shall discuss the 
most robust measurements of this type and point out cer- 
tain advantages that this scheme provides over the usual 
methods of state determination. Realistic experiments 
can be designed along these lines, as we shall see. 

General reasoning. The following counting argument 



already suggests the feasibility of the above idea: The 
number of real parameters to be determined for finding 
p which lives in m-dimensional Hilbert space is rri^ — 1, 
since p is a m x m hermitian matrix with unit trace. On 
the other hand, the most informative repeated measure- 
ments of an observable a) of S are those for which the 
spectrum Ui is non-degenerate; they provide to — 1 inde- 
pendent data, the probabilities tvpfri, i < 1 < m. Thus, 
if measurements are performed on S only, they must deal 
with at least (m? — 1)/(to — 1) = to -f 1 non-commuting 
observables in order to fully determine the unknown p. 
For instance, the 2x2 density matrix p — i(l -I- p- a) of 
a spin-i is parameterized by the to^ — 1 = 3 expectation 

values p — tr( /5 • (7 ) of the Pauli matrices a = (ai , (T2 , (T3) 
with a'i(72 = ifT3, and its determination requires measur- 
ing the spin in m + 1 = 3 non-coplanar directions. Exper- 
iments with three different apparatuses are thus needed. 
This conclusion holds even if we are given the informa- 
tion that the unknown state p is pure (p^ = p or = 1), 
since then the measurement of two components of p does 
not fix the sign of the third one. For to > 3 the construc- 
tion of the needed to -I- 1 non-commuting observables is 
less straightforward, see the discussion in Ref>^. 

In order to design a scheme where p will be determined 
by measurements of a single observable (or, equivalently, 
by commuting measurements only) , it is natural to intro- 
duce * an auxiliary system A that we term the assistant, 
and which lies in a known state f. Let n be the number 
of dimensions of the Hilbert space of A. The compound 
system S-l-A has a mn x mn density matrix TZ = p ® f. 
Ref. ^ proposes measuring one of its observables, a "uni- 
versal quantum observable", VL = 'Y^VLaPa, where the 
spectrum Via is not degenerate so that the eigenprojec- 
tions Pa (with 1 < a < mn, noting that X)q = 1) 
constitute a complete set of mn — 1 commuting observ- 
ables. Such repeated measurements provide mn — 1 inde- 
pendent data Pa — triZPa , the probabilities of the eigen- 
values ^la of Cl. The set Pa also represents the diagonal 
elements of the density matrix p (g) r in the basis that di- 
agonalizes f2. A linear mapping p ^ Pa = ^a.ijPji 
is thus generated, leading from the to^ — 1 real parame- 
ters of p to the mn — 1 data Pa . If the assistant A has 
the same dimensionality to as S, this mapping is repre- 
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sented by a square matrix. In general, if the measured 
observable fl intertwines sufficiently S and A, the deter- 
minant A = det Ra^ij of this matrix may be expected to 
be non-zero. The inverse mapping then solves our prob- 
lem: measurements of Cl, performed repeatedly with a 
single apparatus, yield the probabilities Pa, the knowl- 
edge of which is equivalent to that of p. If n > m, the 
set Pa can still determine /5, but it is overabundant. 

However, the above scheme is not easy to implement 
in practice, since it implies measuring an observable Cl 
(or a commuting set Pa) which thoroughly mixes S and 
A. We propose here a modified procedure, which will al- 
low a much simpler choice for 51. The measurement of 
il is performed not at the time t — 0, at which S is pre- 
pared in the unknown state p and A in the known state 
f, but at a later time t = t. During the lapse < t < r, 
S and A interact, their evolution being generated by a 
known Hamiltonian H . The state of the composite sys- 
tem S+A which is tested is now = UTZoU^, where 
the initial state is TZq = p ® f and the evolution opera- 



tor is U 



-iHt 



The required mixing of p and f being 



thus achieved by dynamics, we can now measure the sim- 
plest possible non-degenerate observable Cl, a factorized 
quantity f2 = cD (g) 6. The observables uj and o of S and 
A have the spectral decompositions to — ^j'^i ^-nd 

6 = X]a=i "aPa ^^'^ projection operator Pa, with 
a = {ia}, takes the form Pa = Pia = t^i ® Pa- Repeated 
measurements of f2, that is, of uj and 6 simultaneously, 
determine the joint probabilities 



Pa = Pia = tr[/(/5® f)C/t(7ri (»p„) 



(1) 



to observe uji for S and Oq for A. (The numbers Pia are 
the diagonal elements of l}{p® f)U^ in the factorized 
basis which diagonalizes a) and 6.) Like above, the map- 
ping p i-^ Pa is expected to be invertible for n > m, 
provided H couples S and A sufhciently. Then simulta- 
neous measurements of lu on the system S and of 6 on the 
assistant A, based on the counting of the events {ia} and 
of their correlations, fully determine p through inversion 
of the equation iQJ . 

For given observables uj and 6 and for a given initial 
state f of the assistant, the precision of this scheme of 
measurement of p relies on the ratio between the exper- 
imental uncertainty about the set Pa and the resulting 
uncertainty on p, which can be characterized by the de- 
terminant A of the transformation For A = it 
would be impossible to determine p by means of Pa- 
The Hamiltonian H and the duration r of the interac- 
tion should thus be chosen so as to maximize |A|. 

Two by two density matrix. In the following we illus- 
trate the above ideas by studying a two-level system S 
(m = 2). We exhibit in particular the best measurement 
schemes that correspond to the largest |A|. We use the 
spin-i representation p = ^{1 + p ■ a). The determina- 
tion of the unknown polarization vector p relies on the 
coupling of S with the assistant A, which we first take 
as another two- level system (n = 2). The observables uj 



and 6 to be measured are the z-components as and S3 of 
S and A, which may be equal to 1 or —1. The projection 



operators are ni 



(1 ± 0-3) and Pa = =•= for i 



and a equal to ±1. Experiments determine the four joint 

probabilities Pa = {P++, P-\ ,P-+,P , } for 0-3 and S3 

to equal 1 or —1. These probabilities are related to the 
three real parameters poip through Eq. , which reads 



Pa 



(2) 
,(3) 



with a — {ia} — {±±} and matrix elements taken in the 
standard representation of the Pauli matrices a and s. 

By construction, the mapping El is such the the 
probabilities Pa are non-negative and normalized for any 
p such that p^ < I. These properties are expressed by 



Ua > \Va\, 



^Ua^l, ^Va=0. 



(4) 



The determinant A of the transformation p 1— s- is 
four times the volume of the parallelepiped having any 
three of the four vectors Va as its sides, e.g., A = 

4t/++ • {v-+ xv^ ). Provided the evolution operator U is 

such that the vectors Va are not coplanar, the transfor- 
mation ((2Jl can be inverted, and p is deduced from the set 
Pa of classical probabilities. Alternatively, p is deduced 
from the knowledge of the expectation values (173), (S3) 
and (17333) at the time t = t, which are simultaneously 
measurable and are in one-to-one correspondence with 
the set Pa. 

We first look for the upper bound of |A| impHed by 
the conditions l@J. First we note that |A| increases with 
\va\ for each a. We therefore maximize A^ under the 
constraints Itfaj = 1 and J2a — 0, that we account 
for by means of Lagrange multipliers A and p. Varying 

iA^ + XY,a Wa\ + pY.a WC find 



A(i/+_ X = A 



(5) 



and other equations resulting from all permutations of 
the a's. This yields symmetric solutions for which the 
four vectors Va, ot = {±±} form a regular tetrahedron: 



Va ■ Vf3 



(6) 



As could have been anticipated, the solutions are not 
unique: they follow one from another by rotation in the 
space of the spins a and permutations of the indices a. 
The corresponding upper bound for the modulus of the 
determinant is |A| = l/(12\/3). 

Let us show that this upper bound can be reached, pro- 
vided the initial state of the assistant is pure. We have 
to construct a unitary operator U that satisfies Eqs. 
EJ, and to find a Hamiltonian H and an interaction time 



T such that U 



-iHt 



Let us exhibit an example of 
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such a solution. We assume that the assistant is ini- 
tially polarized in the z-direction, f — = i(l + (T3), 
and we orient the tetrahedron Va in the direction vj^± = 
(±l,l,±l)/4v^, w_± = (±l,-l,Tl)/4V3. The corre- 
spondence 10) then takes an especially simple form: 



Pi 



= V3( 



P2 



= V3( 



0-3 ) 



P3 



0-3 S3) 



(7) 



yielding directly the density matrix p in terms of the 
expectation values and the correlation of the commuting 
observables (T3 and S3 in the final state. It is easy to 
verify that this correspondence can be achieved under 
the action of the Hamiltonian 



H = ai 



si cos ( 



S3 sm ( 



(S2 - Si 1 sm( 



S3 cos ( 



V2 



,(8) 



where 2<j> — 0.95531 is the angle between v^^ and the 
z-axis, that is, cos 20 = l/VS- Noting that = sin^ 
where x = 1-11069 satisfies cosx = ^cos(f>, we obtain 

U—e^^^'^ = cos(Tsinx) — iiJ sin(Tsinx)/sinx. Taking 
as duration of the evolution t ~ xl sinXj we obtain U = 
cos X — iH. Insertion in (jSJ allows to check Eq. © and 
to get the expected optimal correspondence JT)!. 
The simpler Hamiltonian 



H ^ ^= (Tisi -I- i(s2 sin(^ -I- S3), 



(9) 



which results from |(HJ) through rotation of s, also achieves 
an optimal mapping p ^ Pa, provided S3 is replaced 
by Sismcj) + S3 cos both in the measured projections 
= i(l ± S3) and in the initial state f = The first 
term of Eq. (^J describes, in the spin language, an Ising 
coupling, while the second term represents a transverse 
magnetic field acting on A only. 

Larger assistant. We have optimized above the deter- 
mination of p by coupling S to an assistant A that starts 
in a pure state and has the same dimension n = m = 2 as 
S. It is natural to wonder whether the quality of the mea- 
surement, as expressed by the magnitude of A, may be 
improved for n > m = 2 and/or for mixed initial states of 
A. We may, for instance, consider an assistant consisting 
of q spins, in which case m = 2 and n = 2'. We now de- 
note as S3 some two- valued observable of A which is sub- 
jected to measurement at the time r. The only changes 
in 1^ are the dimension n of the matrix f and the fact 
that the two projection operators Pa no longer constitute 
a complete set in the Hilbert space of A. Experiment still 
provides the four probabilities Pa = tr7?,r(7ri ®Pa) with 
unit sum, where TZr = U{p (8) r)C/^^ is the final density 
operator in the mn-dimensional space of S-t-A, and the 
mapping p ^ Pa keeps the form . The conditions Q 
still hold, since they express simply that the correspon- 
dence or (O preserves the positivity and the normal- 
ization. When obtaining the upper bound 1/(12^3) for 
|A| we relied only on these conditions. Therefore, us- 
ing a larger assistant cannot improve upon the optimal 
solutions found for n = 2 and pure f. 



In all cases, the maximum of |A| is reached for map- 
pings (PlEIl which involve the regular tetrahedron In 
such mappings there exist 4 pure states, pa = ^(1 ~ 4z/a • 
a), for which one probability, Pg, vanishes. 

Completely disordered assistant. Returning to the case 
n = TO = 2 we saw that we could reach the upper bound 
l/(12\/3) of |A| if the assistant is initially in a pure state. 
In order to explore how much is lost if it is in a mixed 
state, we consider the extreme situation in which f ^ ^1 
is the completely disordered state. There is an advantage 
in using such a state, as it is easier to prepare than a pure 
state: one lets the assistant interact with a hot thermal 
bath; for a spin, one leaves it unpolarized. Eqs. ((21 
taken for r ~ show that too simple evolutions may 
lead to a vanishing determinant A. For instance, the 
above evolution which we optimized for f = p+, is 
completely ineffective for f = ^1, since it maps any p 
onto the trivial set of probabilities Pa — j- 

In order to maximize | A| with respect to U for f = ^1, 
we now have to evaluate the vectors Va from We 
only sketch this calculation here. We represent Uai.bj 
as [/+.,+. = VK, [/_.,+ . = WK', [/+._. = VK'X, 
U-..-. = —WKX, in terms of two hermitian positive 
2x2 matrices K and K' such that K"^ + K'^ = 1, and 
of three unitary 2x2 matrices V, W, X in the space 
i — ±, j — ± of S. Using the invariances of the problem, 
we may without restrictions parametrize K — K{9, ip) = 
cos 9 cos (fi + sinOsimpx ■ a with unit vector x, implying 
K' = K{e - I, ip) and - K'"^ = K{2e, 2ip). We char- 
acterize V by V'^a-iV = rf- a, W hy W^^a^W = C ' <5', and 
take X = ^ ■ a, where also ff, C, ^ are unit vectors, with 
X • f = 0. We then obtain Ua — j and the Va, whence 



A = — sin 46 sin 4(p (x x 



{(■Ov + iv-OC -(10) 



We now determine the maximum value of |A| attainable 
for f = il and the corresponding optimal evolution U. 
The best angles are 9 = ip = tt/8. Choosing x = (1, 0, 0) 
and (0, 1, 0), we see that the largest |A| is 1/32. This 
still corresponds to Ha forming a regular tetrahedron, but 
with \va\ = instead of \ija\ = j- Accordingly, for 

the same orientation of the Va's, the factor \/3 in {Tj) is 
replaced by 2. Rather surprisingly, the efficiency of the 
new scheme is not much worse than when the assistant is 
prepared in a pure state. However, the Hamiltonians 
needed now to maximize A require a more complicated 
coupling than in Among various possibilities, we 

present here an example: 



H 



1 



V2, 



1 



CT2S2 



V2 



((Tl + Si) 



(11) 



(the signs of the three terms can be changed inde- 
pendently). This Hamiltonian involves an anisotropic 
Heisenberg interaction and an external field acting sym- 
metrically on S and A. Using Eq. (j2Jl with f = ^1, one 
can check that for r — j{2k + l)n the evolution operator 
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U = c '^"^ leads to an optimal solution with A = 1/32. 
For r = 7r/4 Eq. ((T)) is replaced by 

pi= 2(17353), p2 = 2((s3)cos7+ ((73)sin7), (12) 

P3 = 2(((T3)cos7- (S3)sin7), 7 = ^lli^t.^. (13) 

Conclusion. The non-commutative information con- 
tained in the density matrix of quantum system can be 
transformed by a one-to-one correspondence into ordi- 
nary information associated with a set Pa of ordinary 
probabilities for exclusive events. The price to be paid is 
the introduction of our extra assistant system A. This 
correspondence can be experimentally implemented by 
repeatedly letting S and A suitably interact, then by 
performing each time simultaneous measurements of two 
commuting observables lo and o pertaining to S and A, 
respectively. Counting of events yields the probabilities 
Pia for a) to take the value uji and for 6 to take the value 
Oq. Provided lo is non degenerate so that its number of 
distinct eigenvalues tOi (1 < i < m) equals the dimension 
m of the Hilbert space of S, and provided the number of 
eigenvalues of 6 is also m (at least), the correspondence 
P ^ Pa can be inverted, as we displayed on several exam- 
ples. The second condition implies that the assistant has 
a dimension n > m. Hence an initially unknown p can 
be determined via Pa by means of a single apparatus. 

As compared to the standard determination schemes 
of p based on direct non-commutative measurements of 
S, the present method has several advantages, i) It is 
more economical, since it involves only one observable Co 
of S and one observable 6 of A, whereas direct determi- 
nations require measuring at least m-\-\ non-commuting 
observables of S. ii) This full set of m -I- 1 observables 
is not always accessible in practice. For instance, for a 
two-level atom prepared in some unknown state p, p^ 
is readily measured through the occupation probability 
of the excited state, but pi and p2 can be determined 
only indirectly. Interaction of S with another initially 
known two-level atom A (whose own preparation may be 
straightforward, as we saw above) may provide the full p 
through mere simultaneous measurements of the occupa- 
tion probabilities for the levels of S and A. Hi) It has been 
stressediS that the use of standard statistical and infor- 
mation theoretical methods for dealing with incomplete 
or noisy experimental data cannot be directly extended 
to quantum mechanics, because results are produced 
there by means of non-commutative measurements; in- 
deed, these data pertain to different contexts as they are 
produced by different apparatuses. The present scheme, 
involving only commutative measurements, circumvents 
this difficulty. 

Taking as a criterion of quality of our measurement 
schemes the size of the determinant of the mapping 
/5 I— > Pq, we have explored for m — 2 the conditions 
that lead to the best determination of p for some uncer- 
tainty on the set Pa. For an assistant with dimension 
2, its known initial state f should be pure and the pa- 



rameters of the Hamiltonian should be suitably chosen. 
An example of an optimal evolution is generated by the 
Hamiltonian JH)) or 0. We have also seen that the deter- 
minant cannot be enlarged by use of an assistant with 
dimension larger than 2. (However, 6 may then take 
more than two eigenvalues Oa and the probabilities Pia 
become more numerous than needed; the fact that they 
are related to one another independently of p opens the 
possibility of improving the determination of p through 
cross-check of the data Pa.) 

For TO = 2, the optimal mappings 10, © amount to 
identify, via a dynamical process, the joint probabilities 
Pa for (T3 and S3 with the expectation values in the state 
p of the observables 5 ( 1 — „ ) , where Via = ^ ( 1 ~ • cr) 
pertains to the system S. The 4 observables Cla are pro- 
jection operators, satisfying trfJaii^ = ^ for (3, and 
spanning the space of observables cD of S. For m > 2 
we conjecture that a bound on A may be found by con- 
sidering in the Hilbert space of S a set of projec- 
tions r^Q, satisfying trSla — 1, tv^a^p — 1/(to -I- 1) 
for a 7^ /3, X)q ~ ml, and constituting a ba- 
sis for the observables (i. Then the mapping matrix 
i? in Pq, = Ra,ijPji is expected to be given by 
to(to — l)Ra,ij = Sij — ^a,ij. This form makes the de- 
terminant A stationary under the constraints imposed 
by positivity and normalization alone. As above, there 
are pure states pa — f^a for which one probability, 
Pq, vanishes. This conjecture yields for A^ the upper 
bound to~"'^[to(to -I- 1)(to — 1)^]^~™ , which generalizes 
the TO = 2-result A^ = l/(12\/3)2. 

It might have been expected that a completely disor- 
dered assistant makes the determination of p from ob- 
servation of Pa unprecise or even infeasible; for m = 2 it 
turns out that the best Jacobian in this case is smaller 
than the maximum one only by a factor (\/3/2)'^ ~ 0.65, 
and that Hamiltonians as simple as (|ll|l can be used. 

Indeed, the various types of two-level systems on which 
experiments are currently performed (NMR, quantum 
and atomic optics, spintronics) feature Hamiltonians sim- 
ilar to Q and Hll|l that optimize the process, with Ising 
or Heisenberg types of couplings. For instance, the spin- 
spin interaction between two single-electron quantum 
dots is usually anisotropic due to spin-orbit coupling or 
to a lack of symmetry of the host material; sec^^ for a 
recent discussion. Experiments can therefore easily be 
designed along the above ideas. They will demonstrate 
that the principle of complementarity, which seems to 
imply that different measurement devices are needed to 
fully determine a quantum state, can be by-passed by 
using an assistant, even completely disordered. 
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